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Abstract 

The dynamics of self-gravitating fluids is analyzed within the framework of a collisionless Boltzmann equation 
in the presence of gravitational fields and Poisson equation. Two cases are analyzed: a system with baryonic and 
dark matter in a static universe and a single system in an expanding universe. The amplitudes of the perturbed 
distribution functions are considered as a linear combination of the collision invariants of the Boltzmann equation. 
For the system of baryonic and dark matter, the Jeans mass of the combined system is smaller than the one of 
the single system indicating that a smaller mass is needed to initiate the collapse. For the single system in an 
expanding universe it is not necessary to make use of Jeans ’’swindle”and it shown that for small wavelengths 
the density contrast oscillates while for large wavelengths it grows with time and the Jeans instability emerges. 


1 Introduction 

The search for structure formation from gas clouds is an old subject that goes back to 1902 when Jeans [1] used 
the system of phenomenological equations of mass and momentum densities together with the Poisson equation 
and showed that small perturbations in the mass density, pressure, velocity and gravitational potential in a static 
background with wavenumber smaller than the Jeans wavenumber could evolve with time. In terms of balance 
of forces the fluctuations can grow in time if the inwards directed gravitational force is larger than the outwards 
directed internal pressure of the gas. 

Jeans theory describes the gravitational instability of self-gravitating systems by searching for conditions that 
small perturbations can grow and leads to a collapse of the system [2, 3, 4, 5, 6, 8, 9, 10]. It was formulated 
before the knowledge of the Universe expansion and one has to take into account the Jeans “swindle”, which 
imposes that the Poisson equation is valid only for the perturbations, since the background solution of constant 
mass density, pressure, gravitational potential and vanishing velocity satisfy the balance equations of mass and 
momentum densities, but not the Poisson equation. 

The analysis of the small perturbations by using the same phenomenological equations of Jeans but taking into 
account the expansion of the Universe was due to Bonnor [7] in 1957. The background solutions in this case satisfy 
the balance equations of mass and momentum densities and Poisson equation. In this context there is no necessity 
to invoke Jeans “swindle” (see e.g. [2, 4, 7]). 

Another way to examine Jeans instability is to consider the system of equations composed by the Newtonian 
version of the Boltzmann equation and Poisson equation (see e.g. [4, 5, 11, 12, 13, 14]). In these works it is 
considered a static Universe and Jeans “swindle”. The dispersion relation follows by taking into account small 
perturbations of the equilibrium values of the distribution function and gravitational potential. 

Today it is well known that the matter content of the Universe is composed by baryonic matter which comprises 
atoms of all categories and dark matter which does not emit or interact with electromagnetic radiation. The ratio 
baryonic matter to dark matter estimated is about one to five. Dark matter plays a key role in structure formation 
because it feels only the force of gravity: the gravitational Jeans instability which allows compact structures to 
form is not opposed by any force, such as radiation pressure. As a result, dark matter begins to collapse into a 
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complex network of dark matter halos well before ordinary matter, which is impeded by pressure forces. Without 
dark matter, the epoch of galaxy formation would occur substantially later in the universe than is observed. 

The aim of this work is to obtain the collapse criterion from the Newtonian version of the Boltzmann equation 
and Poisson equation for a system of baryonic matter and dark matter and for a single system in an expanding 
spatially flat Universe. In both cases the amplitudes of the perturbed distribution functions are taken as a linear 
combination of the collision invariants of the Boltzmann equation. For the system of baryonic and dark matter it is 
shown that the mass needed to initiate the collapse is smaller than the one of a single constituent and that the ratio 
of the dispersion velocities have a significant role in the structure formation. For the single system the equilibrium 
distribution function is written in a comoving frame which takes into account Hubble’s law. Furthermore, the mass 
density is a solution the Friedmann and acceleration equations for a pressureless fluid. There is no necessity to use 
Jeans “swindle”, since the equilibrium distribution function and gravitational potential satisfy both the Boltzmann 
and Poisson equations. For this case the Jeans instability is connected with the growth of the density contrast for 
large wavelengths, which is a parameter related with the local increase of matter density. 

2 Basic equations 

In the phase space spanned by the space and velocity coordinates (r, v) a state of a non-relativistic gas is char¬ 
acterized by the one-particle distribution function / = /(r,v,t), such that f{r,'v,t)(Prd^v gives the number of 
the particles in the volume element (Pr about the position r and with velocity in the range Pv about v. The 
Boltzmann equation governs the space-time evolution of the one-particle distribution function / in the phase space. 
In the presence of a gravitational potential $ and in the absence of collisions between the particles, the Boltzmann 
equation is given by (see e.g. [16, 15]) 


SJ + v-V/-V$-9v/ = 0. (1) 

Here the force per unit mass which acts on a particle is only of gravitational nature F = —Vd>. 

The Boltzmann equation (1) follows also from a relativistic version of it in the presence of gravitational fields, 
where a one-particle distribution function / = /(r,p,t) in the phase space spanned by the space coordinates and 
momenta (r, p) satishes the Boltzmann equation (see e.g. [17]) 

p^dj-r^,p^p’^dp.f = o. ( 2 ) 

Above the mass-shell condition p^p^ = rrP(? - with m denoting the particle rest mass - was considered. In the 
non-relativistic Newtonian limiting case —>■ me, p —>■ mv and the Christoffel symbol Fqq —>■ so that (2) 

reduces to (1). 

The connection of the gravitational potential with the one-particle distribution function is given by Poisson 
equation 


V2(j) = AttGp = JttG J mfPv, (3) 

where G denotes the gravitational constant and p is the mass density of the fluid. 

The analysis of cosmological problems is based on the cosmological principle which states that the Universe 
at large scales is spatially homogeneous and isotropic. The metric which satisfies this principle is the Friedmann- 
Lamaitre-Robertson-Walker metric. For a spatially flat Universe the line element in this metric is written as 
= (edt)^ — a(t)^ (dx^ + dy'^ -\- dz^) , where a{t) is the cosmic scale factor. 

Einstein’s field equations for a Universe dominated by a perfect fluid reduce to two coupled differential equations 
(see e.g. [2, 3, 4]) 
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SttG 


3 


P, 



(4) 


which are known as the Friedmann and acceleration equations, respectively. Above the dot denotes the derivative 
with respect to time and are the mass density and the pressure of the sources of the gravitational field, 
respectively. 
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In this work we are interested in matter dominated Universe, where the sources of the gravitational field are 
pressureless p = 0. In this case equations (4) can be solved for the mass density as a function of the cosmic scale 
factor, yielding 


P = Po 



(5) 


where po and ag are the values of the mass density and cosmic scale factor at t = 0, respectively. 

Another relationship which will be used below is Hubble’s law r = (d/a)r which relates the recession velocity 
of an object r with its physical distance r in an expanding Universe. The quantity = a/a is known as Hubble 
parameter. 


3 Jeans instability for systems of one component 


We start by analyzing the Jeans instability for systems with one component. This subject was discussed in some 
books, see e.g. [4, 5], but here a new methodology is introduced, which is based on the collision invariants of the 
Boltzmann equation [15]. The collisionless Boltzmann (1) and Poisson (3) equations imply a coupled system of 
equations for the distribution function /(r,v,t) and gravitational potential <I’(r,t), which we reproduce below: 

dtf + V • V/ — V4> • dvf = 0, = 47 rG / mfcfv. ( 6 ) 


At equilibrium the distribution function and the gravitational potential of a self-gravitating system do not depend 
on time and the collisionless Boltzmann equation ( 6 )i is satisfied if the distribution function depends only of the 
velocity of the particles /o(v) and the gravitational potential gradient vanishes V<I>o =0. Note that the condition 
V'I’o = 0 may follow from symmetry considerations, because in a homogeneous system there is no preference in the 
direction of the gravitational potential gradient. The same restriction for the gravitational potential gradient is also 
considered for a system described by the balance equations of mass density and fluid velocity, since constant values 
of pressure and density and vanishing fluid velocity and gravitational potential gradient satisfy these equations. 
However, the condition V$o = 0 does not satisfy the Poisson equation (6)2 due to the fact that its right-hand side 
refers to the mass density of the self-gravitating system. In order to remove this inconsistency one makes use of 
the Jeans “swindle”, which considers that the Poisson equation is valid only for the perturbed distribution function 
and perturbed gravitational potential. 

We require that the equilibrium state is subjected to small perturbations /i(r,v,t) and $i(r,<): 

= fo{v)[l + h{r,v,t)], $(r,t) = $ 0 -k $i(r,t). (7) 


The equilibrium value of the gravitational potential is considered as a constant in order to fulfill the condition 
V'I’o = 0 and the equilibrium distribution function is the Maxwellian one 


/o(v) 


m (271(7^)3/^ ’ 


( 8 ) 


where a = y^kT/m is the dispersion (thermal) velocity and k, T denote the Boltzmann constant and the temperature 
of the system, respectively. Note that the collisionless Boltzmann equation ( 6)1 is a differential equation for the 
distribution function and its solution depends on the problem we are analyzing. Here we are interested in analyzing 
disturbances which can be represented by plane waves so that we have required that the distribution function is a 
sum of an equilibrium distribution function plus an arbitrary function which characterizes a plane wave solution of 
small amplitude. 

We insert the representations (7) into the Boltzmann and Poisson equations ( 6 ) and get the following system of 
equations for the perturbed quantities h and $ 1 : 


/o [9t/r + V • V/r] - V$i • Sv/o = 0, 


= AttG 


mfohd^v. 


(9) 


In ( 9)1 we have neglected products of V<I>i with h and since h and $1 are supposed as small quantities. 
Furthermore, we have considered that the Poisson equation refers only to the perturbed values of the distribution 
function and gravitational potential and taken V4>o =0. 
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The perturbations h and $i are now represented by plane waves of frequency uj and wavenumber vector k such 
that 


/i(r, V, t) =/ii(v) exp [i (k • r — wt)], <I)i(r, t) = (jiexp [i (k • r — wt)]. (10) 

In the above equations ^ is a constant amplitude while hi(y) is given as a combination of the collision invariants 
of the Boltzmann equation, namely: 1, v and (see e.g. [15, 16]). Its expression is written as 

/ii(v) = A + B • V + Hv^, (II) 

where A, B and D are unknown constants. 

The insertion of (10) and (II) together with (8) into the perturbed Boltzmann and Poisson equations (9) lead 
to 


(A + B • V + Hv^) [w — k • v] — k • V— = 0, (j) + AttG {A + Sa^D) p = 0, 


( 12 ) 


where k = jkj. 

Now a system of equations ioi A, B = H ■ k, D and </> can be obtained from the multiplication of (12)i by the 
collision invariants (l,v,v^) and integration of the resulting equations, yielding 


oj [A + 3a^D) — a^B = 0, ojB — 


A + 5a^D+^ 


= 0, UJ {3A + 15a'^D) - = 0. 


(13) 


The system of equations (13) together with ( 12)2 has a non-trivial solution if the determinant of the coefficients 
A, B, D and </> vanishes, which implies the dispersion relation: 


w* (1 — k^ + uj'l') — 0. 

Above the dimensionless wavenumber modulus k^, and frequency w* are given by the relationships 


k* = 


v'IttG'p k,7 


— 1—) w* — 




(14) 


(15) 


Here we have introduced the adiabatic sound speed Vg = \/5/3ct and the Jeans wavenumber kj = ^JAnGpIvs- 
We have the following solutions for the dimensionless frequency from the dispersion relation (14): w* = 0 and 


w* = ± 



(16) 


From the above equation we may infer that for k > kj (small wavelengths), w* is a real quantity and the perturbations 
- due to the factor exp {—iujt) - will propagate as harmonic waves in time. However, if k < kj (big wavelengths), 
w* is a pure imaginary quantity and one of the perturbations will grow with time, while the other will decay with 
time. The perturbation which grows with time is associated with Jeans instability. 

A very simple model is described in the literature to understand the Jeans instability. Consider a mass M 
enclosed in a spherical volume of radius A in which there exists a mass density inhomogeneity. This inhomogeneity 
will grow if the gravity force Fq per unit mass is greater than the opposed pressure force Fp per unit of mass, i.e.. 


Fg 


GM 


A2 


oc 


GpX^ 

A2 


GpA > rp OQ —oc —, 
A 


since 


v^g oc 


p 

p' 


(17) 


Now identifying the Jeans wavelength in terms of the Jeans wavenumber Aj = 27r/kj = g jAnGp and A = 27r/k, 

the inequality k < kj follows. 
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4 Jeans instability for systems with two components 


In this section we shall analyze Jeans instability by taking into account two collisionless Boltzmann equations, one 
for the baryonic matter and another for the dark matter which are connected with the Poisson equation. We shall 
use the indices b and d for the baryonic and dark matter, respectively. The corresponding evolution equations for 
the distribution functions of baryonic matter /{, = /(x, Vh,t) and dark matter fd = f{'K,Vd,t) are given by the 
collisionless Boltzmann equations 


dtfb + Vb • V/b - • d^Jb = 0, dtfd + Vd • Vfd - v$ • d^Jd = 0, 


which are connected with the Poisson equation 


= JttG mbfbd^Vb + J rridfdd^Vd^ = 47rG(ph + pd). 


(18) 


(19) 


Above {mb^rnd) and {pb,Pd) are the masses and mass densities of the baryonic and dark matter, respectively. 

Following the same methodology of the previous section we suppose that the distribution functions /(x, Vf,,t), 
f{x,Vd,t) and the potential $ are subjected to small perturbations hb{v,Vb,t), hdir,'Vd,t) and <l>i(r,t) from their 
equilibrium values /°(vb), /^(v^) and $ 0 , respectively. Hence we write 

/(i’,v&,t) =/°(vb) [1 + /it,(r,Vb,t)], f{r,Vd,t) = f^{vd)[l + hdir,vd,t)], $(r, t) = $0 + $i(r, t). (20) 


where the equilibrium distribution functions are the Maxwellians 


/°(GV&,t) 


Pb e 

mb (27r(T^)3/2 ’ 


/d(i'Wd,0 


Pd e 

md (271(7^)3/2 


( 21 ) 


Here ab = a/ kTb/mb and ad = ^JkTd/md are the dispersion velocities of the baryonic and dark matter which are 
associated with their temperatures Tb and Td^ respectively. 

From the insertion of the representations (20) into the Boltzmann (18) and Poisson (19) equations yields 

/° [dthb + vf, • Vhb] - V$i • = 0, /° [dthd + vd • Vhd] - V$i • d^,fd = 0, (22) 

= JttG mbfbhbd^Vb + J mdf%dd^v^ . (23) 


Here we followed the same methodology which was used to obtain the system of equations (9) and considered the 
Jeans ’’swindle”. 

We represent the perturbations as plane waves of frequency w and wavenumber vector k 

hb{r,Vb,t) = hl{vb)exp (ik • r - iujt) , hdir,Vd,t) = /i;)(vd)exp (ik • r - iut) , $i(r,t) = (/)exp (ik • r - iujt) , (24) 

where the amplitude <j) is constant, while hi and hi are given in terms of the collision invariants of the Boltzmann 
equations (l,Vb,v^) and (l,Vd,v^): 

f^bi'^b) = Af, + Bf, • Vb + Db'vf, kdi'Vd) = Ad + B(i • Vd + (25) 

Above Ab, Ad,Jib,T^d, dDb, Dd are constants 

Now the insertion of (24) together with (25) and (21) into (22) and (23) leads to 

(Aq, + Ba • Va + UaV^) [w - k-Va] - k • = 0, a = b,d (26) 

k^cj) + AttG [{Ab + ialDb) Pb + {Ad + ia^Dd) Pd] = 0. (27) 

In order to get a system of algebraic equations for A?,, A^, = Bf, • k, Bd = ^d- k, Db, Dd and (j) we multiply (26) 

by the collision invariants (l, v?,, v^) and (l, v^, v^), respectively, and integrate the resulting equations, yielding for 
(a = b, d) 


uj (A„ + 3alDa) - Baal = 0, 


ujB„ — 


Aa + 5a Da H—y 

CTaJ 


= 0 , 


w (3A, + 15alDa) - 5Baal = 0. (28) 
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A non-trivial solution of the system of equations (27) and (28) is found if the determinant of the coefficients 
Ah, Ad, Bh, Bd, Db, Dd, 4> vanishes, which implies the dispersion relation: 




^2 

, .2 , 

H - 9 

[k^ 1 

L Pd 

V \ 

^d 



= 0 . 


(29) 


Here we have introduced the dimensionless wavenumber k* and frequency w* defined in terms of the dark matter 
Jeans wavenumber kj = ^JAirCpd/v'^ - with rif = ^JhJiud denoting the dark matter sound speed - since as it was 
explained in the introduction the dark matter begins to collapse into a complex network of dark matter halos well 
before ordinary matter. The dimensionless wavenumber and frequency read 


k kuf w 

kj i/47rGpd ’ * y/AirCpd 


(30) 


We note that the dispersion relation (29) is a function of two ratios pd/pb and Od/ob- The mass density ratio 
Pd/Pb can be associated with the present value of the density parameter ratio Pd/Pb = ^d/^b ~ 5.5 [18], since this 
ratio has not changed considerably during the evolution of the universe. For the dispersion velocities ratio Od/cfb 
there is no fixed value. One value we shall use is taken from Ref. [19] where Milky Way-like galaxy simulations 
including both baryonic and dark matter were performed. From this work we have inferred that in one of the 
simulations where Maxwellian distributions are considered the ratio is given by (Jd/(Jb = 170/93 « 1.83. 

From the dispersion relation (29) we have the following solutions for the dimensionless frequency: 



Without the baryonic matter, i.e. for pb = ab = 0, the dispersion relation (29) reduces to: 

ujt+[l- k^] 00^ = 0, (33) 

which has the solutions w* = 0 and w* = i-^k^/kj — 1, i. e., we recover the Jeans solution for one component. 

Two solutions of (31) provide real values for w* for any k, so that we have harmonic waves in time for these two 
solutions. However the other two solutions furnish imaginary values for w* for some values of k* so that instead of 
harmonic waves the amplitude of the disturbance will grow or decay. The one which grows is associated with Jeans 
instability. The value of k* where w* changes from the imaginary value to the real value is obtained by taking w* 
equal to zero in (29) and we get 


k* 




Af 


(34) 


This equation is interpreted as the ratio of two Jeans wavenumbers, the one kj*^ refers to the system dark-baryonic 
matter while kj to the dark matter. 

We can also analyze the amount of mass which is necessary to initiate the collapse, which is the Jeans mass 
contained in a sphere of diameter equal to the wavelength A = 27r/k. Hence, we can build the ratio of Jeans masses, 
one for the system dark-baryonic matter Mj^ and another for the dark matter alone Mj, namely 


Pd / 



(35) 


In Table 4 the ratio of the Jeans masses of the system baryonic-dark matter and dark matter are given as 
functions of the dispersion velocities ratio for fixed values of the mass densities ratio. From this table we infer that 
if we increase the dispersion velocities ratio, the mass needed to begin the collapse becomes smaller in comparison 
with the mass where only one constituent is present. This can be understood, since when the ratio Ud/ub is large 
the dispersion velocity of the baryonic matter is smaller than the one of the dark matter and the baryonic matter 
hardly overcome the escape velocity of a given gravitational field. 
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1.00 

1.20 

1.40 

1.60 

1.83 

2.00 

2.20 


0.9199 

0.8338 

0.7481 

0.6662 

0.5791 

0.5206 

0.4585 


Table 1: Ratio of Jeans masses as function of the ratio of the dispersion velocities Udjoh for pdjPh = 5-5. 


5 Jeans instability in an expanding Universe 


In this section we shall analyze Jeans instability by taking into account the collisionless Boltzmann and Poisson 
equations (9) but in an expanding Universe where a pressureless fluid is the source of the gravitational field. 

Now the equilibrium distribution function is written in a comoving frame as 


/o(r,v,f) 


P 1 
m (27r(T2)3/^ 




(36) 


thanks to Hubble’s law f = {ala)v. Note that according to (5) the mass density p is only a function of time. 

The gravitational potential function 

O-rr 

=—Gpr\ (37) 

and the equilibrium distribution function (36) satisfy the Poisson and Boltzmann equations 


V^$o = 47rG 


mfocfv 


AirGp, 


dtfo + V • V/o - V$o • Sv/o = 0, 


(38) 


by taking into account the Friedmann and acceleration equations (4) and provided that the dispersion velocity 
is proportional to the inverse of the cosmic scale factor and does not depend on the spatial coordinates, namely 
a{t) oc l/a{t). 

As in the previous sections we require that the equilibrium state defined by the distribution function (36) and 
gravitational potential (37) is subjected to small perturbations characterized by h{r,-v,t) and $i(r,t) such that 

/(uv,t) =/o(r,v,t) [1 +/i(r,v,t)], $(r, t) = $ 0 ( 1 -, 0 + 4’i(r, t). (39) 

Note that (39) is similar to (7), the difference between them lies in the dependence of the equilibrium distribution 
function and gravitational potential at equilibrium, which now are functions of the space-time coordinates. 

Now the insertion of (39) into the collisionless Boltzmann and Poisson equations (6) leads the system of equations 
for h and $ 1 : 


/o \dth -\- V • V/i — V<i>o ■ d^h] — V$i • 9v/o = 0 , = JttG J mfohcfv. (40) 

Note that the underlined term above does not show up in (9), since it was supposed that V$o =0. Above the 
products of V$i with h and dvh were also neglected. 

The perturbations h and $1 represented by plane waves of wavenumber vector q/a(t), reads 

q- 

a{t) j ’ 

The factor l/a(t) in the wavenumber reflects the fact that the wavelength is stretched out in an expanding Universe. 
Here we cannot assume a harmonic wave in time, since the factors of the equations depend on time. Hence, (j)(t) 
is an amplitude that depends on time, while hi{r,v,t) is given as a combination of the collision invariants of the 
Boltzmann equation in a comoving frame, i.e., 1, v — and (v — ^r) : 

(42) 

where A{t)^ functions of time. 


/ii(r, v,t) = A{t) + B{t) 


a 

V-r 

a 


D{t) 


<I>i(r, t) = exp i 


q- r 

l{t) 


(41) 


/i(r, V, t) = hi{r, V, t) exp i 
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The methodology used here is the same as in the previous sections. We insert (41) and (42) together with (36) 
and (37) into the perturbed Boltzmann and Poisson equations (40) and multiply the resulting equation from the 
Boltzmann equation by the collision invariants (l) v ~ fr, (v — and integrate. Hence we get the following 


system of differential equations: 


dA{t) , o 2dDit) 

„dA{t) 


dt 


15(7 


dt 
2d D{t) 
dt 


+ i — B{t) — 6-a^D{t) = 0, 


dB{t) .q 


2 r 


dt 


A{t) + 5a'^D{t) 


4>{t) 

/t2 


--H(t)=0, (43) 


■i5—B{t) - 30-a^D{t) = 0, 


;(j){t) + AnG [A{t) + 3a'^D{t) p = 0. 


(44) 


Above we have introduced B{t) = B(t) • q. 

In galaxy formation is usual to introduce the density contrast which is a parameter that indicates where there 
are local increase in the matter density. It is defined by the ratio of perturbed and unperturbed mass densities 
6p = 'p! p- Here the density contrast reads 5p = A{t) + 3a‘^D{t), and we can obtain from the system of differential 
equations (43) and (44) that A is a constant (say A = 1) and the following differential equation for the density 
contrast 



rd’p 



4A5 

25AQr§ 


= 0 , 


where 


27rao 
Aq —-, 

q 


IOttcto 

3^/4:^TGpo 


(45) 


In the above equation we have used the relationship a = ao(67rGpo^^)^^^ which follows from the integration of (4). 
Furthermore, t = ty/GTrGpo is a dimensionless time and the primes are derivatives with respect to r. The solution 
of (45) is given in terms of Bessel functions of first kind (A/rs) and reads 


dp = T 6 


Gi Js 



+ G2«/_5 
2 




1 + 


5t3 


where 



(46) 


The solution (46) differs from the one that comes from the analysis of the balance equations for mass, momentum 
and entropy densities in an expanding universe (see e.g. [4]) due to the underlined term which refers to the time 
evolution of the density contrast in a matter dominated universe. For large values of the ratio Aj/Aq (small 
wavelengths) the Bessel functions imply oscillations of the density contrast while for small values of Aj/Aq (large 
wavelengths) the first term in (46) furnishes two solutions one decaying according to 1/t and another growing with 
However due to the underlined term with the density contrast will grow with time for small values of 
Aj/Aq and this solution corresponds to Jeans instability. In Fig. 5 it is shown the time evolution of the density 
contrast for different values of A, i.e. for different values of the ratio Aj/Aq. 


6 Conclusions 

In this work we have analyzed Jeans instability using the collisionless Boltzmann and Poisson equations for a system 
of baryonic and dark matter in a static universe and a system of a single constituent in an expanding universe. 
For the system baryonic-dark matter the Jeans mass of the system is smaller than that of a single component 
indicating that the structures are formed earlier than in the latter case. Furthermore, the dispersion velocities 
ratio dark-matter/baryonic-matter has influence of the Jeans masses, since small dispersion velocity of the baryonic 
matter in comparison with the one of the dark matter implies that the baryonic matter hardly overcome the escape 
velocity of a given gravitational field. For the single component in an expanding Universe the background solutions 
for the distribution function and gravitational potential satisfy both Boltzmann and Poisson equations without 
the introduction of Jeans “swindle”. The equilibrium distribution function is written in a comoving frame that 
takes into account Hubble’s law and the mass density is a solution of Friedmann and acceleration equations for 
a pressureless fluid. The Jeans instability is connected with a parameter that indicates the local increase in the 
matter density, the so-called density contrast. It is shown that for large wavelengths the density contrast grows 
with time. 
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Figure 1: (color online) Density contrast 5p as function of time r for different values of the ratio A = -^/IS/SAj/Aq. 
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